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Abstract 

Event-by-event fluctuations of hadronic patterns in heavy-ion collisions are studied 
in search for signatures of quark-hadron phase transition. Attention is focused on a 
narrow strip in the azimuthal angle with small Ay. The fluctuations in the gaps be- 
tween particles are quantified by simple measures. A scaling exponent a is shown to 
exist around Tc- An index ^ is shown to characterize the critical fluctuation; it is a 
numerical constant ^ = 0.05 it 0.01. All the measures considered in this gap analysis 
are experimentally observable. Whether or not the theoretical predictions, based on 
simulations using 2-dimensional Ising model, are realistic for heavy-ion collisions, anal- 
ysis of the experimental data suggested here should be carried out, since the existence 
of a scaling behavior is of interest in its own right. 
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1 Introduction 



If a quark-gluon system created in a heavy-ion collision at high energy undergoes a second- 
order phase transition to hadrons, one should expect large fluctuations in the hadron densities 
from region to region. Thus in any small interval of hadronization time the particles produced 
on the surface of the system, where the temperature T is at the critical temperature Tc, 
should form patterns of clusters and voids of all sizes. The question is how such fluctuations 
can be detected. Since the detector integrates over the hadronization time, the hadron 
patterns formed at different times, when added, tend toward a uniform sum. In Refs. Q 
we have studied the problem and found a way to overcome the difficulty. On the one hand, 
one must make severe cuts in the transverse momentum pj- (accepting particles only in a 
narrow window of width Ap^), while on the other hand, appropriate measures that are 
insensitive to the randomization process following hadronization should be used to quantify 
the critical fluctuations. We proposed a void analysis that identifies the scaling properties 
of the fiuctuations. The suggested measures in the analysis are in the two-dimensional space 
of rapidity (y) and azimuthal angle (0). The analysis has been applied to the NA49 data 
of Pb-Pb collisions at CERN-SPS 0, but remains to be considered for the data collected at 
BNL-RHIC. 

In this paper we aim at simplifying the experimental problem of detecting the critical 
fluctuation, while making one of the theoretical parameters in the void analysis more exper- 
imentally relevant. They are accomplished by reducing the 2D space to ID so that a void 
becomes a gap. The study of gaps in rapidity was undertaken by us earlier in connection 
with multiparticle production in hadronic collisions [Q. The proposed gap analysis has been 
applied to the NA22 data and provides a way of testing the predictions of the dynamical 
models on soft processes. We now apply the gap analysis to the phase transition problem 
in heavy-ion collisions. Besides showing what can be expected theoretically, we believe that 
the suggestions made here are more amenable to the experimental situation. An important 
outcome of this study is the identiflcation of a numerical constant that characterizes the 
critical fluctuation. It can be checked by experiment. 

Our theoretical tool, as in |T|, 0, is still the simulation of critical phenomenon on the 
2D Ising lattice. Not only is the simulation simple to execute, it is also physically relevant, 
since the quark-hadron phase transition is very likely to be in the same universality class as 
the Ising system 0. We map the 2D Ising lattice to the surface of the quark-gluon plasma 
cylinder formed in a central collisions of heavy ions. We select a row in the lattice and map 
that to a strip on the plasma cylinder conflned to a narrow interval in rapidity. Our proposed 
analysis is for the experimental examination of the fluctuation of gaps between particles in 
the azimuthal variable. 

2 Gap Analysis on the Lattice and for the Experiments 

The relationship between the Ising problem and the hadron density problem has been dis- 
cussed in detail in Refs. Q. We summarize here only the essentials for what is needed in 
this paper. 

We start with a 2D lattice of size L^, where L = 288. On the lattice we deflne cells of 
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size e^, where e = 4. The density pi at the zth cell is defined to be 

A = Ac2^(q) , (1) 
where Cj is the net spin at cell i defined to be positive along the overall magnetization, i.e., 

Ci = sgn^(TJ^(Tj, (2) 

\jeL2 / jei 

aj being the lattice spin at site j, and A an unspecified factor relating to the particle 
density in y-cj) space. Near Tc, pi fiuctuates from cell to cell; however, the existence of 
correlation of all length scales leads to clusters of all sizes formed out of neighboring cells 
having non- vanishing p,. 

In Refs. [|I], 3 we considered bins of size 5^, with 5 being integral multiples of e. We 
regarded a bin to be empty if the average bin density, p^, is less than a threshold density po- 
We then defined a void to be a collection of contiguous empty bins, connected by at least 
one side between neighboring empty bins. The measure used to quantify the fiuctuation of 
void sizes was shown to be insensitive to the value of po. 

Without undercutting the importance of doing the void analysis on experimental data 
with high statistics, especially at high energies, we want now to consider in this paper the 
possibility of doing a similar analysis in ID, since if clustering of hadrons occurs in 2D, they 
should also be visible in ID. Not only will the statistics improve, there will also be no need 
to do binning to define voids, as we shall see. 

Let us consider a row of C cells, where we initially take C = L/e = 72. We map that 
row to a slice of the plasma cylinder in the central region with width Ay in rapidity. What 
Ay is will be discussed below. We also will limit pt to a narrow interval Apy, as discussed 
extensively in Refs. 0. Thus we have only the azimuthal angle as our space of variable. 
The row of C cells is mapped to the interval < < 27r. 

We simulate N^. configurations (where A^^g = 3 x 10^) on the Ising lattice in 2D at T = Tc, 
and determine the density pi at each cell in each configuration. The simulation cannot 
be done in ID because a ID Ising system does not exhibit critical behavior. From a 2D 
configuration we choose an arbitrary row and examine the values of pi in each cell with 
1 < i < C . The value of A in Eq. (0) is effectively set to 1 in what we do with pj, since we 
always compare it to a threshold po in units of A. We define an occupied cell to be one in 
which 

A > Po , (3) 

and we place one (and only one) particle at cell i. If pi < po, then it is an empty cell. In 
mapping to the variable, we place a particle at (pi randomly in the interval (z — l)27r/C < 
< z27r/C, when the zth cell satisfies Eq. (^. Thus after the whole row is mapped to the (f) 
variable, there are particles in the space, where fiuctuates from event to event. An 
event corresponds to a configuration on the lattice. 

On the lattice the average number of occupied cells (A^) in a row, after averaging over 
all configurations, depends on the threshold po. In the experiment the average number of 
particles (A^) in the space depends on AyApx- In our analysis we shall vary po. The 
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predicted dependence of our measure on po is to be checked by experiment by varying 
AyApT- In that sense the variable po is not a theoretical parameter devoid of experimental 
relevance. 

In an experiment of heavy-ion collisions the single-particle inclusive distribution in (p, 
dn/d(f), may not be uniform, because the impact parameter is not always exactly zero (even 
if a centrality cut is made), or because the detector efficiency is not symmetric in 0, mostly 
likely both. It is then better to use the cumulative variable X [|^, ^, defined in our case by 

^w=r^#7ri^#'. (4) 



The range < < 27r is thus mapped to < X < 1, and the density of particles in X, 
dn/dX, is uniform. The distribution dn/d(f) is determined after many events, but the variable 
X is used for exclusive distributions event-by-event. In a simulation dn/d(f) is uniform, if N^. 
is large; we nevertheless use the X variable between and 1 just so that our proposed gap 
analysis corresponds to what is to be done with experimental data. 

Consider an event with N particles in our window, labeled by i = 1, ■ ■ ■ , X, located in 
the X space at Xj, ordered in accordance to Xj < Xj+i. Now define the distance between 
neighboring particles by 

Xi = Xj+i -Xi, i = 0, ■ ■ ■ , X, (5) 

with Xo = and X^r+i = 1 being the boundaries of the X space. We call these Xj's gaps, 
which clearly satisfy 

E^. = l- (6) 

1=0 



Let us define for each event the moments 

1 



N 



tE^I (7) 



which, for positive q, puts more emphasis on the large gaps than the small ones. The set 
{Gq} with 2 < q < Q, Q being some number less than 10, say, will be our quantification of 
the event structure. The set fluctuates from event to event, especially at the critical point. 
Our method is to use {Gq} as a basis to construct a measure that can reveal the critical 
behavior. 

Among the various erraticity measures considered in Ref. we choose 5*^ for our con- 
sideration here, since it has the simplest behavior. Starting with the definition 

Sq = {Gq\nGq), (8) 

where the average (■ ■ ■) is performed over all events, 5*^ is defined by 

S, = Sq/Sf (9) 

where s^* is as defined in Eq. (||) but with the statistical contribution to Gq only. Thus the 
deviation of Sq from 1 is a measure of how strong the dynamical fiuctuations are relative 
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to the statistical fluctuations. In our simulation we generate random numbers between 
and 1 to populate the X space for the calculation of G^*. In an experiment one can do 
something similar to generate a random event. 

We shall study the dependence of Sq on q for various values of po, using configurations 
simulated on the Ising lattice. The same can be done with the experimental data from heavy- 
ion collisions. If hadrons are formed in the latter due to a second-order phase transition, 
then the experimental Sg should reveal some features that are similar to our theoretical Sg 
to be presented below. 

3 Results on Critical Fluctuations 

Using A'^e = 3 X 10^ configurations simulated on the Ising lattice at T = Tc = 2.315 in units 
of J/ks where J is the near neighbor coupling and fc^ the Boltzmann constant, we have 
determined Sq for a range of po- For po = 20 the dependence of Sq on q is shown in a log- log 
plot in Fig. 1. The solid line is a straight-line fit of the calculated result, giving strong 
evidence for the power-law behavior 

Sq OC q''. (10) 

The value of the exponent is a = 1.67±0.02. The fact that In Sq deviates unambiguously from 
implies that Sq is a statistically significant measure of nontrivial dynamical fluctuation. The 
power-law behavior is not necessarily a consequence of the dynamics of critical phenomenon, 
since similar behavior has been found before in the case of hadronic collisions . However, 
the exponent a is indicative of critical fluctuations; it is an order of magnitude larger than 
the value of a = 0.156 obtained for pp collisions at y/s = 20 GeV. 

To have one exponent a for all g is a very economical description of the critical fluctu- 
ations. We can then investigate the dependence of a on po. Before so doing, let us first 
find a replacement of po by some quantity that is directly measurable. In the preceding 
section we have related po to the intervals AyApT as the variable under experimental con- 
trol that can be used to tune the average multiplicity (A^) accepted in the (p window. Thus 
(A^) is a quantity that is both experimentally observable and theoretically computable on 
the lattice. It is therefore a suitable replacement for po. Nevertheless, we prefer to use an 
even better one that is the average number of gaps (M), which is also observable. On the 
lattice in accordance to our convention of counting gaps in Eqs. and (^, we have simply 
M = A^ + 1. However, it should be recalled that in our simulation we have adopted the 
rule that, when p, exceeds po, only one particle is placed in a cell, not more, no matter how 
high Pi is. That procedure makes sense in view of our measure Gq of the event structure, 
since a tightly packed cell with many particles in it would have very small gaps that make 
negligible contribution in Eq. (0). Experimentally, if there are particles whose momenta are 
indistinguishable, or nearly so, whether they are separately counted or not also makes no 
significant difference in the calculation of Gq. Thus to allow for such possibilities it is better 
to count the number of distinguishable gaps, rather than the number of particles. 

On the lattice we have simulated the configurations for a range of po from 20 to 200. For 
every value of po we can calculate the average number of gaps, (M). Figure 2 shows the 
dependence of (M) on po at Tc. The straight-line segments are just the linear interpolations 
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between neighboring points of po where the calculation has been made. With such a definitive 
relationship at hand, any quantity that depends on po will in the following be shown as a 
function of (M), so as to render it amenable to experimental verification. 

For each of the higher values of po examined, we have found power-law behavior of Sq, 
as in Fig. 1. Thus the exponent a can be determined in each case. In Fig. 3 we show the 
dependence of a on (M). Remarkably, the dependence is very linear. If we parameterize it 
as 

a = ao + e(M), (11) 

we obtain 

ao = -0.258, ^ = 0.055. (12) 

The nature of critical fiuctuation is now seen to be reduced to a simple formula, Eq. ([TT|) , 
when the moments of gaps are used to describe the event structure. If we further put em- 
phasis on the property that is independent of (M) , then the slope ^ in Eq. (0) emerges as a 
numerical output of the theory that relies on no numerical input. This is perhaps the most 
succinct characterization of the critical phenomenon, beside the critical exponents. The lat- 
ter depend on the temperature of a critical system near Tc. In heavy-ion collisions T is not 
directly observable, so the corresponding critical exponents cannot be measured experimen- 
tally (see, however, Ref. [^). Here we have an index ^ that is eminently measurable and is 
the only numerical constant that can be meaningfully associated with critical fiuctuation. 

So far our study has been done only at T = T^. To see how the results change when T 
deviates from T^, we repeat the above analysis for a range of T. Figure 4 shows Sq vs q in 
log-log plot for 2.27 < T < 2.80 in units of J/ks- Note that linearity is quickly lost when T 
goes below = 2.315. For T > the linearity persists for a limited range 2.315 < T < 3.1, 
but the slope is reduced. For T > 3.2 the dependence bends over at high q and the linearity 
is lost. For the range of T where a can be determined from linear fits, we show in Fig. 5 
a(T) up to T = 2.8. The sharp peak that occurs at T = Tc is a remarkable manifestation 
of the critical behavior. For T less than Tc, so many hadrons are produced that the gap 
distribution rapidly tends toward the statistical. For T larger than Tc, fewer hadrons are 
produced, and it takes more T — Tc difference for the statistical fiuctuation to dominate. If 
T were experimentally controllable, the measurement of a(T) would be an excellent way to 
determine the critical temperature of the quark-gluon system. That being not the case in 
reality, we can only learn from Figs. 4 and 5 that, if the system hadronizes in a range of 
T around Tc, the most significant portion of the contribution to Sq would come from the 
immediate neighborhood of T = Tc, and that only the a value around Tc is experimentally 
relevant. 

The above analysis is done for po = 20. We can, of course, repeat the analysis for other 
values of po. In each case we find the peak value of a at Tc. At each T where a can be 
meaningfully determined, we can investigate its dependence on (M), just as we have done 
in Fig. 3 at T^. In each case a linearity is found that allows us to determine the slope index 
^ in Eq. (|Tl|). The result is shown in Fig. 6, where ^(T) also exhibits a peak at Tc. We can 
now conclude that if hadronization is to occur around Tc, we expect the measurable value of 
^ to be around 0.05. 
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Having established the properties of the observables as functions of T due to the dynamics 
of the critical system, we now investigate the question of stability with respect to changes in 
the size of the detector window, which is kinematical. We have so far considered a lattice of 
size 72x72 cells, from which we choose one row of C = 72 cells from each configuration. We 
now want to vary the length and width of the row. We shall do so by setting T at Tc. First, 
we consider a row of C = 54 cells from each configuration, and later a row of C = 36 cells. 
In each case it does not matter whether the row is mapped to a correspondingly shorter 
range in since dn/d(j) is converted to dn/dX in the X space before the gap analysis is 
performed. However, the average number of gaps (M) does change, so the simulation with 
shorter rows does correspond to a shorter ranges of in the experiment. In Fig. 7 we show 
the result on a vs. (M) for the C = 72, 54, and 36 cells. The straight lines are linear fits, 
which are evidently very good. The values of the slope ^ are : 

^ = 0.055 ± 0.005, C = 72 

0.057 ±0.005, C = 54 

0.048 ± 0.007, C = 36 (13) 

The case of C = 36 cells yields a slightly lower value and relatively larger errors, as is 
reasonable, since there are fewer points in a shorter range. Nevertheless, all three values of 
^ agree within errors. We therefore conclude that ^ is stable against variations in the length 
of the window. 

Next, we consider variations in the width of the window. In our simulation so far we 
have used only rows with 1-cell width. In an experiment the width Ay of a window in 
will have to be adjusted in order to vary (M). To check the stability of our result against 
variations in the width of our strip on the Ising lattice, we consider combinations of two and 
three rows. Let k be the row index, and i the cell index in the row, as before. Thus the 
location of a cell on the lattice is now denoted by ki. We take r adjacent rows and combine 
them by performing vertical average of r cells at each horizontal position i , i. e., we define 

1 

A = - I] Pki (14) 
k=i 

and then proceed in the gap analysis using pi as in the preceding section. It is important that 
the r rows be adjacent so that the short-range vertical correlation of the cells can influence 
the value of so as to reflect its dynamical content; otherwise, if the r rows are randomly 
chosen, the averaging in Eq. tend to render pi more statistical and its distribution 
in i more uniform. In short, taking r adjacent rows corresponds to widening Ay in the 
experiment. In our calculation we consider only the case r = 2 and 3. Note that we take 
the average in Eq. (p!4D and then vary po to change (M). We could have added pki (without 
dividing by r) and not vary po! that is equivalent to averaging p^i and dividing po by r. Our 
chosen procedure allows more continous change in po. 

The result of our analysis for r = 1, 2 and 3, and C = 72, are shown in Fig. 8. Roughly 
speaking, the slopes are essentially the same within errors, so ^ may be regarded as inde- 
pendent of r. More precisely, we find 

C = 0.055 ±0.005, r = l. 
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^ = 0.044 ±0.003, r = 2, 

^ = 0.045 ±0.003, r = 3. 



(15) 



The slight decrease in ^, as r increases, is not unreasonable, since averaging over r rows 
tends to suppress dynamical fluctuations. In the limit r becoming very large, clearly only 
statistical fluctuations remain. We may summarize Eqs. (pTS]), ( plSj) and Fig. 6 by the grand 
average 

^= 0.05 ±0.01. (16) 

If this can be verified by experiments, then one should be able to claim that a signature of 
critical transition has been observed. If the value of ^ is not confirmed, yet the power-law 
behavior of Eq. (|T0|) is shown to exist in the data, with or without the linear dependence 
of a on (M) in Eq. ([TT|) , that would still be an exciting experimental finding, suggestive of 
dynamical fiuctuations. 

There is the usual question about final-state interaction and the dilution of the dynamical 
signature due to randomization. The answer depends on the type of measure for that signa- 
ture. The issue has been addressed in Ref. 0, where the dependence on the number of steps 
of final-state scatterings is examined, in Ref. |]l| where configuration mixing is considered, 
and in Ref. where different options in simulating time evolution are investigated. What 
one learns from all those studies is that the measures considered are not strongly affected by 
the final-state randomization. Since the gap analysis is a derivative of the void analysis [0, 0, 
the same conclusion follows here. The one reminder that we should emphasize is that the 
window ApT in the transverse momentum should be kept as small as possible to minimize 
the overlap of particles emitted at different times. 

4 Conclusion 

We have studied the problem of even-to-event fluctuations of the hadronic patterns in phase 
space in heavy-ion collisions in search for detectable signatures of second-order quark-hadron 
phase transition. We have reduced the complication of voids in two dimensions to the simpler 
problem of gaps in one dimension. Using the moments of gaps to construct an entropy-like 
measure Sq, we have found a power-law dependence on q with an exponent a. It is the 
dependence of a on the average number of gaps that yields the index ^, which serves to 
characterize critical fluctuation. We have found the stability of ^ = 0.05 ± 0.01 against 
variations in the length and width of the detector window in (p. When the temperature of 
the system is moved away from Tc, the power-law behavior of 5*^ on q persists in a narrow 
range of T around Tc, and the values of a and ^ show strong peaks at Tc. Thus is a measure 
of the critical behavior and is a number that arises out of the study of fluctuations without 
any numerical input. It is highly signiflcant that C, can be checked by experiments, since all 
measures leading to its determination are designed to be observable. 

In heavy-ion experiments it is not difficult to make cuts in Apx and Ay to limit the 
average number of particles in a narrow strip in (p to the range from 10 to 40. If those 
particles are found to be non-uniformly distributed for every event, the gap analysis proposed 
here is a way to quantify those fluctuations. When the exponent a and the index ^ are found 
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to exist, there are numerous variables under the control of the experiments to vary. We can 
mention, for example, the position of the y slice, the value of p^, the centrality of collisions, 
the nuclei sizes and the cm. energy. It would be very interesting to see the dependence 
of ^ on the total transverse energy Et, since it can provide us with some idea of when the 
critical behavior is lost as the nuclei overlap gets to be too small to create a quark-gluon 
plasma. Also, for nonzero impact parameter, even if critical fluctuations exist, the index ^ 
may depend on which sector of the (j) space the analysis is performed, since there is no 
invariance in elliptic flow. One can envision a rich variety of phenomenological studies once 
the exponent a and index ^ are found in the experimental data. They can provide valuable 
information about the quark-gluon system. The application of the gap analysis to the data 
is therefore strongly urged. 
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Figure Captions 



Fig. 


1 


Sq vs q in log-log plot for T = and po = 20. 


Fig. 


2 


The average number of gaps (M) vs po at Tp. 


Fig. 


3 


The exponent a vs (M) at Tc. 


Fig. 


4 


Sg vs g in log-log plot for a range of T at po = 20. 


Fig. 


5 


q; vs T at Po = 20. 


Fig. 


6 


^ vs T at Po = 20. 


Fig. 


7 


a vs (M) at Tc for various number of cells in a row. 


Fig. 


8 


a vs (M) at Tc after averaging over various number of rows. 
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